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METHOD TO DETERMINE PROPERTIES OF A SEDIMENTARY BODY 
FROM THICKNESS AND GRAIN SIZE DISTRIBUTION AT A POINT 

WITHIN THE BODY 



FIELD OF THE INVF.TVTTmvr 

[0001] This invention relates generaUy to the field of geophysical prospecting. 
Specifically, the invention is a method for detennining flie size, shape, and grain size 
distribution of a sedimentary body fi-om a measurement of its thickness and grain size 
distribution at one point 

BACKGRQTTNn OF THF TIWFTVTTOivr 

[0002] In the oil and gas industry, seismic prospecting techniques commonly are 
used to aid in the search for and evaluation of subterranean hydrocarbon deposits. A 
seismic prospecting operation consists of three separate stages: data acquisition, data 
processing, and data interpretation, and success of the operation depends on 
satisfactory completion of all three stages. 

10003] In the data acquisition stage, a seismic source is used to generate an 
acoustic signal that propagates into the earth and is at least partially reflected by 
subsurface seismic reflectors. The reflected signals are detected and recorded by an 
array of seismic receivers located at or near the surface of the earth, in an overlying 
body of water, or at known depths in boreholes. 

[0004] During the data processing stage, the recorded seismic signals are refined 
and enhanced using a variety of procedures that depend on the nature of the geologic 
structure being investigated and on &e characteristics of the raw data. In general, the 
purpose of the data processing stage is to produce an image of the subsurface from the 
recorded seismic data for use during the data interpretation stage. 
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[0005] The puipose of the data interpretation stage is to detennine infonnation 
about the subsurface geology of the earth fiom the processed seismic data. The 
results of the data inteipretation stage may be used to detennine the general geologic 
structure of a subsurface region, or to locate potential hydrocarbon reservoirs, or to 
guide the development of an already discovered reservoir. 

I0006J At present, the conclusions that can be made after the data inteipretation 
stage are generaUy limited to broad descriptions of the size, shape, and general nature 
of subsurface structures. These descriptions may, for example, provide an indication 
of the total volume of hydrocarbons, which might be retained in such structures. 
However, present technology does not allow the analyst to detennine preferred 
locations within a stmcture for drilUng weUs, except in a very broad sense. In 
addition, when an exploration well has been driUed, present technology does not 
allow an analyst to be able to accurately characterize the nature of the subsurface 
stmcture in locations other than the immediate region surrounding the well. 

[0007] As wiU be understood from this summary of the data interpretation stage 
of the seismic prospecting operation, it is desirable to be able to predict the broad 
nature of a subsurface stmcture of interest using only seismic data and a minimal 
number of exploration wells, preferably just one exploration well. Such a capabiUty 
would faciUtate estimation of hydrocarbon volume in place and production rates early 
in the hydrocarbon exploration and development process. As will be understood to 
those skilled in the art, the hydrocarbon volume and rate of production depend on a 
variety of factors, including the grain size distribution of the sand deposit that makes 
up the reservoir in which the hydrocarbons are found. 

SUMMARY OF THE TNVENTTON 

[0008J A method to detennine the properties of a water-lain sediment body from a 
measurement of grain size distribution and deposit thickness at one location in the 
body is disclosed. In one embodiment, the method comprises (a) detennining the 
flow properties at the measured location, (b) extrapolating the flow properties back to 
the mlet through which the depositing flow was emitted, (c) detennining at least one 
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property of the water-lain sediment throughout the sediment body by modeling the 
flow properties using the extrapolated flow properties at the inlet from step (b) as a 
boundary condition. The flow properties associated with the sediment body include 
flow velocity, suspended sediment volume fractions, deposition time, and flow height. 
The properties of the water-lain sediments include, in addition to the flow properties 
associated with deposition of the sediments, the thickness of the sediment body, the 
size of the body, the shape of the body, and the grain size distribution at each point 
within the body. 

[0009] A second embodiment of the method to determine the properties of a 
water-lain sediment body from a measurement of grain size distribution and deposit 
thickness at one location in the body is also disclosed. The method comprises (a) 
estimating the flow height at the measured location, (b) determining the flow 
properties at the measured location from the grain size distribution and deposit 
thickness at the measured location, (c) extrapolating the flow properties and deposit 
properties along a characteristic curve intersecting the measurement point, (d) 
identifying the maximum exfr^olated deposit thickness, (e) identifying the inlet point 
along the characteristic curve as the point where the flow properties are consistent 
with the inlet conditions, (f) repeating steps (a) through (e) until the maximum 
extrapolated deposit thickness is consistent with the flow properties at the identified 
inlet point, (g) predicting deposit properties throughout tiie whole sedimentary body 
by using the extrapolated flow properties at the identified inlet point as boundary 
conditions for a forward model. 

BRIEF DES CRIPTION OF TTTF. ntt SWINGS 

[0010] Figure 1 is a plan view of a fluid flow which is depositing a sedimentary 
body, including the flow boundaries; 

[0011] Figure 2 is an elevation view corresponding to plan view Figure 1; 

[0012] Figure 3 is an elevation view corresponding to plan view Figure 1 after 
deposition has occurred; 
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f0013J Figure 4 is a flow chart of an embodiment of the invention; 
(00141 Figure 5 is a flow chart of a second embodiment of the invention; 

DETAILED PESCRinPTTOM 

(0015] hi the foUowmg detailed description, the invention will be described in 
comiection with its preferred embodiment. However, to the extent that the foUowing 
description is specific to aparticular embodiment or a particular use of the invention, 
this is intended to be illustrative only. Accordingly, the invention is not limited to the 
specific embodiment described below, but rather, the mvention mcludes all 
alternatives, modifications, and equivalents falling within the true scope of the 
appended claims. 

mi€l The present invention is a method for predicting size, shape, and internal 
properties of a sedhnent deposit fiom a weU penetration of the deposit at one location 
without havmg to drill multiple wells. SpecificaUy, the inventive method can hi 
appUed to finding the grain size distribution at every point within a water-Iam 
sediment deposit as well as the shape and size of the deposit using a measurement at 
one pomt of the deposit thickness and grain size distribution. 

[0017] As will be miderstood to those skilled in the art, the deposition of clastic 
sedimentary bodies typically begins with a flow of sediment-laden water from a 
confined channel, such as a river mouth, into an open region, such as a basm The 
pomt where the sedhnent-laden flow enters an open region where deposition occm^ is 
known as the inlet. Initially such flows expand freely and deposit sediment as the 
flow decelerates. Thereafter, as the deposited sediment grows in height, the deposited 
sediment begms to obstruct the flow field. EventuaUy. the deposit becomes 
sufficiently large that the flow is diverted around liie deposit. This results in a new 
path and mlet for the flow field to an open region beyond or adjacent to the old 
dqx>sit nie deposition process Ihen repeats, and a second body m the system is 
created, hi addition, more than one such body may be actively built within the system 
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at a time. OveraU, the process produces a deposit consisting of stacks of sedimentary 
bodies, which is a useful model for the structure of hydrocarbon reservoirs. 

[00181 In the present inventive method, appUcants recognized that since each 
body stops building when it is large enough to divert the flow, the scale of a body is 
related to the flow that builds it. Other properties of the body including its size, 
shape, and the internal properties, such as grain size distribution, are also controlled 
by the flow. 

[0019] As will be understood by those skilled in the art, the flow characteristics 
which create such bodies can be analyzed &om Navier-Stokes and continuity-based 
fluid flow principles appKed to a sediment laden flow beneath stationary clear water 
that can be entrained, such as turbidity currents in the deep ocean. Such analyses can 
be extended to flows without vertical entraimnent, such as a shallow water layer under 
air. and these extensions are also within the scope of the present invention. 

[0020] Applicants recognized that applying reasonable assumptions to these 
principles facilitates an analysis procedure, which can be tied to a physical and 
empirical principle-based quantification of the characteristics of fluid flow-deposited 
sand bodies. This tie between fluid flow principles and sedimentary body deposit 
characteristics faciUtates a method, which allows characterization of the overaU 
properties of the deposit based on only a minimum amount of initially available data. 
That method, in its various embodiments, is the subject of the present invention. 

[0021] The inventive method utiHzes an assumption that the grains which 
compose the deposit are carried and deposited by a sediment-laden fluid flow with 
height hCx,y), x-y^locity u(x,y), and j;-velocity vCx,y). The fluid flow contains a 
volume fraction Ci(x,y) of grains of the ith grain size. The grains in the ith grain size 
bin all have a nominal characteristic diameter of rf, and fell through still water with a 
terminal velocity (also referred to as a settling velocity) v^^. Furthermore, the method 
is built on the assumption that the sedhnent, both in the flow and in the deposit, may 
be characterized by bins of a common nominal grain diameter. Each such bm may, 
for example, be characterized by a nominal diameter d,, and mclude grains havmg ari 
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actual diameter which ranges fiom 84% to 1 19% ofdi. The time duration over which 
dq)osition persists to form the deposit is T. and the resulting deposit thickness is 
z(x,y). 

[0022] Figures 1, 2, and 3 depict the assumptions and parameters used in the 
present invention. Figure 1 depicts a plan view of fluid flow 10 with flow boundaries 
12 and 14. Inlet 15 for fluid flow 10 is ceutered, for convenience, at the origin of the 
X and y-axes, and flow emitted from the inlet moves initially in the positive x 
direction. At the inlet 15, the flow boundary has an initial width 17 and half-width 8 
and expands in the positive x direction. Also depicted is an outline of the deposit 16 
formed by the flow, 

[0023] Figure 2 depicts an elevation view corresponding to plan view Figure 1. 
The fluid inside the flow boundaries 12 and 14 of Figure 1 is comprised of two layers. 
Figure 2 illustrates the two layers of fluid as a clear layer 20 above a sediment-laden 
layer 28. The sediment-laden layer 28 is also referred to herein as the turbid water 
layer. The elevation of the bottom topography before the deposition process occurs is 
24. The height of the sediment-laden water layer at inlet 15 is 26. The height of the 
sediment-laden water layer 28 varies based on location, as evidenced by the fact that 
at a different point in Figure 2, Uie sediment-laden water layer has a different height 
27. 

[0024] Figure 3 is an elevation view corresponding to the plan view of Figure 1 
after deposition has occurred. As in Figure 2, the fluid flow in Figure 3 is depicted as 
being comprised of clear layer 20 above a turbid layer 28. The elevation of the 
subsurface topography after deposition is 30. This elevation consists of the elevation 
of the original bottom topography 24 in addition to the thickness of the newly 
deposited sediment layer 36. 

[0025] Table 1 provides a detailed Ust of the parameters used in the present 
method. As is indicated in Table 1, the height ho and velocity uo of the turbid water 
flow at the inlet are assumed to be constant at the inlet, across the width of the inlet 
along the j;-axis. The velocity uo is also constant verticaUy through the thickness of 
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the flow. The inlet half-width A is also defined. Thereafter, based on the equations 
described herein, the height h and velocity components « and v of the turbid water 
layer are functions of position with respect to the x and j; axes. Velocity components 
« and V are assumed to be constant over the depth of the turbid water layer. The turbid 
water layer is assumed to have a volume fraction Q of grains of the jth grain size 
averaged through the depth of the layer, but that volume fraction varies with jc and j; 
throughout the flow. 

[0026] The deposit, which is formed from the flow of the turbid water layer, is 
assumed to be the result of the net deposition from the flow. The net deposition rate 
from the flow of grains of the /th grain size is the difference between the deposition 
rate A of grains of the /th grain size from the flow onto the bottom and the erosion 
rate Et of grains of the /th grain size from the bottom back into the turbid water layer. 
The deposition process varies with the sediment concentration, and the erosion 
process varies with the flow velocity and deposit grain size distribution, so the net 
deposition wiU also vary with ^ and j; locations. The net deposition process is 
assumed to be in steady-state; in other words the flow velocities, heights, and 
sediment volume fractions of the tuibid water layer at each jt and location are 
constant throughout the deposition time T. The deposit which results is assumed to 
have a constant grain size distribution d through the depth of the deposit, but that 
distribution can vary with and location throughout the deposit. The result of the 
deposition process after deposition time J is a deposit whose thickness z varies with x 
imdy. It is assumed that deposition persists for deposition time T, and that the flow is 
uninfluenced by the growing body it deposits until tune T. At time T, it is assumed 
that the body has reached a critical height, and the process of building the body is 
complete. Physically, this critical height is the height sufficient to divert the flow 
away from the body so that it is no longer building the body. 

[00271 A first embodiment will now be described. With reference to Figure 4, 
this embodiment involves procedures to detennine the flow properties at the inlet and 
to use the flow properties at the inlet as boundary conditions to model properties 
throughout a sedimentary body. The properties include the size and shape of the 
sedimentary body, the flow field that deposited the body, and the grain size 
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distribution at any point within the sedimentary body. As illustrated in Figure 4, the 
flow properties at the measurement location are determined (step 401). The flow 
properties determined in step 401 are extrapolated to the inlet location (step 402). 
Finally, the extiiapolated flow properties at tiie inlet are used as boundary conditions 
for determining the internal properties of the body by modeling the flow field which 
built the deposit and the associated deposit properties (step 403). The individual steps 
will be described in greater detail in tiie foUowing paragraphs. 

[0028] First, the flow properties of the sedimentary body are determined. The 
thickness of the deposit and grain size distiibution can be determined by well core 
data, seismic data or modeling of a sedimentary basin. In the present invention, both 
the thickness of the deposit z and the grain size distribution of the deposit (?, are 
assumed to be a function of both tiie deposition A and tiie erosion Ei fliat occurs 
during the fluid flow process. More specifically: 

Dq>osit thickness: 

z{x,y)^ A(x,y,T) - A(x,y,0) = ^j:iD,ix,y)-E,(x,y)) [i] 

Grain Size Distribution in Deposit: 

G, (x,y) = ^ ^>(^.y)-E,(x,y) 

[0029] Now referring to Figure 4, step 402, tiie flow properties are tiien 
extrapolated to tiie inlet location. The modeling of tiie depositing flow can be 
performed if flie boundary conditions are known. An outflow boundary condition is 
applied along tiie edges of tiie simulation region, except for tiie j-axis as seen in 
Figure 1. Along flie j;-axis, a u=0 boundary condition is appHed except at tiie inlet, 
where tiie boundary conditions must be determined as in tiie inventive metiiod. Those 
sMUed in tiie art will recognize that otiier boundary conditions could be selected 
around tiie simulation region excluding the inlet, and otiier choices of non-inlet 
boundary conditions are wifliin tiie scope of tiie invention. 
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[0030] Now refening to Figure 1, step 403, the flow field is modeled using the 
flow properties at the inlet as a boundary condition. In this model of the fluid flow 
process, deposition is assumed to be a function of the following relatively simple 
relationship to the settling velocity vj/ of the grains in still water: 

Deposition Rate: 

where, for example, vsi may be calculated using the settling velocity equations 
pubhshed by Dietrich (1982). The factor yo is the factor by which the sediment 
concentration at the bottom of the flow exceeds C,, the verticaUy averaged sediment 
concentration in the flow. Typically, yo = 2, though other choices of yo are also within 
the scope of the invention. 

[0031] It is generaUy understood that the erosion that occurs during the fluid flow 
process follows a complicated, empirically derived, relationshqj to the velocity of the 
flow. One such empirical erosion relationship is specified by Garcia and Parker, 
1991: 



Erosion Rate: = min 



where, Z, = A— J?^, 



Pso) 



^ 3 

V 

Iq these equations, the term. 



[4b] 



[4c] 



[4d] 



is a shear velocity which is a function of the magnitude of the local flow velocity and 
a drag coefficient/ « 0.001-0.01 depending on the nature of the bottom. Typically,/ 
= 0.002. The grain size distribution is d, and the grain diameter associated with the 
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50th percentae of grain size in the deposit is pso. The paiameter p = 1.3 x 10"' and A, 
= l-0.288<r^ where tr^ is the standard deviation of the grain size distribution measured 
in the logarithmic "phi" units familiar to geologists. 

[0032] The flow field at all points can be computed if the boundary conditions, 
particularly the flow conditions at the inlet, are known. The term "flow field" is used 
to mean the flow velocity parameters u(x,y) and v(x^), as well as the height of the 
turbid water layer h(x,y) and sediment volume firactions Ci(xj>J. The flow equations 
derive &om the Navier-Stokes equations and continuity equations for a sediment 
laden flow beneath stationary clear water, and are important because the shape and 
internal properties of a deposited body depend on the flow field that deposited the 
body. For the motion of a turbid water layer, Parker, et al. (1986), Imran, et al. 
(1998), and Bradford and Katopodes (1999) all derived an appropriate system of 
equations. In a steady-state form, these equations can be written as foUows: 

X-Momentum: 
Y-Momentum: 

Sediment Conservation: 

dhuCi 9hvC , 



Fluid Conservation: 
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dhu dhv r-^ t 

The multipUer tenn on the right hand side of the fluid conservation Equation 8a is, 
g 0.00153 

0.0204 + F^-* ' [8^] 

and will be recognized to those skilled in the art as an eatrainment coefficient often 
given the symbol e^. Persons skUled in the art will recognize other variations of this 
term that are also within the scope of this invention. This tenn is an empirical 
function of the local Froude number at each point in the flow. The Froude number is, 

Froude Number: 



[8c] 



As indicated. Equations 5 through 8a derive from the Navier-Stokes and continuity 
principles for a sediment-laden flow beneath stationary clear water. However, certain 
assumptions are made in using these equations to determine the flow field in the 
present invention. Persons skilled in the art with the benefit of the disclosures herein 
wiU recognize other equations could be used to calculate the flow properties. All such 
equations are intended to be within the scope of this patent. These assumptions 
include: 1) that the sediment-laden portion of the flow is assumed to have at each x,y 
location, a uniform velocity and fixed sediment concentration profile vertically 
through its thickness, and 2) that the flow field is constant across the entire width of 
the inlet, as ejqpressed in the following boundary conditions: 

Inlet boundary conditions: 



u(0,y) = 



"o \y\ ^ b 

0 \y\>b' t9a] 
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[0033] As wiU be understood by those skiUed in the art, fluid flow systems are 
often preferably analyzed using a characteristic curve coordinate system. Such 
analysis begins with the definition of the characteristic curve relationships. 



~^hu -^ = ^v, [10] 

where s is the charactCTistic curve parameter. 
[0034] In this coordinate system. 



X~Momentum: 



& "~2^^~Sb^ «k+^.>)V«'+v^ , [11] 



Y-Momentum: 

dv 1 aQA" f X ri — T 

~d^^~2^~^ v(e,+c^)V«^+v^ . [12] 

[0035] Define w as the velocity magnitude so that. 



dw _l ( du dv^ 

'd^~w["~d^^''-^)- [13] 



Using Equations 11, 12, and 13, 



Velocity Magnitude: 
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dw Rg dCrh^ w V 

^ = — [14] 

The inlet boundary condition from Equation 9a can be expressed as. 
Inlet boundary condition: 



1^0 M^b 

0 \y\>b- [153 



[0036J The sediment mass conservation Equation 7 is transformed into ihe 
characteristic curve system: 

Sediment Mass Conservation: 

dC. / 

It is sometimes useful to write d as the sum of Ceu an equilibrium tenn for which the 
dq)osition and erosion rates exactly cancel, and Cp,, the excess sediment volume 
fraction. 

Equilibrium profile: 

Equation 16 for Q then gives. 
Sediment Excess: 

[0037] The fluid conservation equation along the family of curves indexed by s 
will be. 
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Fluid Conservation 
dh 



ds'^^-"^' " \Qx^^y [19] 

The curves indexed by s are not quite characteristic curves of the fluid conservation 
equation, due to the divergence term that remains in Equation 19. Assume that this 
term can be written as a function of the variables along the characteristic curves. 
From this function, the flow properties may be determined throughout the 
sedimentary body. M this specific case, assume that the function is zero. However, 
persons skilled in the art with the benefit of the disclosures herein will recognize 
alternative functions of the variables along the characteristic curves. AU such 
functions are intended to be within the scope of this patent. 

(du dv\ , 

ta7"'5;J = -^(^'^'^'^') = 0 [20] 

Under that assumption. Equation 19 can be rewritten as: 
Approximate Fluid Conservation: 
dh , 

[21] 

[0038] Using Equations 1, 16, and 21, the velocity magnitude from Equation 14 
becomes, 

Approx. Velocity Magnitude: 

dw_^ri-^z ^ 

ds 2 [^w'^"-^- [22] 

[0039] A second embodiment of the present invention will now be described. 
With reference to Figure 5, this embodiment involves an iterative analysis procedure 
to establish the flow parameters at the inlet and then forward model the flow from the 
inlet to estabUsh the full size, shape, and internal properties of the body. First, the 
flow height at the measurement location is estimated or guessed (step 501). Next, the 
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flow properties (such as flow velocity, suspended sediment volume fiactions, and 
deposition time) are determined at the measurement location from the measured grain 
size distribution and deposit thickness (step S02). The flow height, flow velocity, 
suspended sediment volume fractions, and deposit thickness are extrapolated along 
the characteristic curve (step 503). The maximum extrapolated deposit thickness is 
identified (step 504). The mlet point is identified as the point where the extrapolated 
flow height, flow velocity, and suspended sediment volume fractions are consistent 
with an empirical inlet condition (step 505). The maximum extrapolated deposit 
thickness and the extr^olated flow properties at the identified inlet point are 
compared with a theoretical relationship between maximum deposit thickness and 
inlet flow properties. If the extrapolated values are not consistent with the theoretical 
relationship, then steps 501 - 506 are repeated, adjusting the estimate of flow height in 
step 501 until consistency is achieved (step 506). The extrapolated flow properties at 
the identified inlet location are used as boundary conditions for a forward model to 
predict the deposit shape, size, and grain size distribution at any point (step 507). 
This embodiment will be discussed in more detail in the foUowing paragraphs. 

[0040] The initial step of this embodiment of the present method of Figure 5, step 
501, is to estimate the flow height at the measurement point. The initial estimate 
might typically be one meter. The estimate is iteratively refined by: extrapolating the 
flow field away &om the measurement point, checking consistency conditions of the 
exti^polated results, and updating tiie flow height estimate until the consistency 
conditions are met. 

[0041] In the second step of this embodiment of the present invention. Figure 5, 
step 502, the flow field above the measurement point is determined from the deposit 
tiiickness and grain size distribution at the measurement point. There are many 
combinations of flow field properties that would produce the same grain size 
distaibution and deposit thickness at the measurement point. Additional assumptions 
may be introduced to determine a unique set of flow properties. One such method is 
described here. Those skilled in tiie art will understand that other assumptions could 
be ^pUed to constrain the choice of flow properties. Such other methods for 
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detennining the flow field above the measurement point are also within the scope of 
the invention. 

[0042] Assume that C« is slowly changing so that Equation 18 is equal to zero. 
Solving for Cpt in Equation 18 then gives. 



Sediment Excess: 



[0043] The measured deposit thickness and grain size distribution at the 
measurement location constrain the choice of flow properties. In tenns of Cjn, 
Equation 1 is. 

Deposit thickness: 

T 



■^Z^'LroVs.Cp, , [24] 



and Equation 2 is. 

Grain Size Dist, in Deposit: 

i 

The grain size distribution G/ gives no information about the grain size fi-actions in the 
flow which do not deposit, except that their volume fi:action in the flow is less than or 
equal to the equilibrium volume firaction for that grain size. For this reason, the total 
volume fiaction of non-depositing grains wiU be combined into the single parameter 
Cs, so that. 

Total Sediment volume: 
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Where the summation in Equation 26 is over those grain size bins for which Gi is non- 
zero. 

I0044J The unknown properties of the flow at the measurement point are then C«, 
w, C& and T/(l-4,). The flow height h was fixed in the first step (step 501). These 
parameters are determined by adjusting w, Cs, and 77(7-^> until the resulting C« from 
Equation 23 best satisfies the constraints of Equations 24 and 25. Experiments with 
this embodiment suggest that the depth constraint of Equation 24 be used first to 
define a one-to-one relationship between w and T/(l-^), and then a second constraint 
such as Equation 25 then be used to select a single pair of w and T/fl-^ values from 
the one-to-one relationship. Persons skilled in the art will recognize that the bedding 
type observed in the core fix>m which the grain size distribution is measured could 
alternatively be used to estimate w, providing the second constraint Persons skilled 
in the art will fiffther recognize that other methods could be employed to estimate 
flow velocity from the grain size distribution, providing the second constraint. Such 
other methods for determining flow velocity are within the scope of this invention. 
The value of T/(l-^) thus detemiined is then used along with the measured thickness 
and grain size distribution to calculate €„ in the flow from Equations 24 and 25. The 
value of w thus determined is then used along with the measured grain size 
distribution to determine Cei in the flow from Equation 17. Experiments with this 
embodiment suggest that Cs may typically be ignored and assumed to be zero. 



[00451 In the third step of the embodiment of the method. Figure 5, step 503, the 
flow field and deposit lliickness are extrapolated away from the measurement point. 
Note that Equations 16. 21. and 22 describe the flow along aU characteristic curves. 
The measurement point lies on some characteristic curve. Since the equations are 
identical on all characteristic curves, it is not necessary to determine which curve the 
measurement point is on. The measurement location can be arbitrarily assigned the 
location s=0, and integration proceeds in the negative s direction to move toward the 
inlet along the characteristic curve and in the positive s direction to move away from 
the inlet along the characteristic curve. Persons skilled in the art wiU recognize that 
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the extrapolation along the characteristic curve could also be accomplished by other 
solution methods, including, but not limited to. control volume methods [Patankar 
1980] for solving simultaneously the flow properties at aU points along the 
characteristic curve. 

[0046] In the fourth step of the embodiment of the method. Figure 5, step 504, the 
maximum deposit thickness obtained along the characteristic curve is identified. 

[0047] In the fifth step of the embodiment of the method. Figure 5, step 505, the 
point along the characteristic curve is found where the extrapolated flow field is 
consistent with a specified inlet condition. The inlet condition is a criterion that 
aUows the location of the inlet to be determined fiom the extrapolated flow properties 
along the characteristic curve. For example one suitable inlet condition would be, 

ho-aW, [27] 
The corresponding mlet width can also be calculated according to the relationship, 

b = a^u^ 

0 [28] 
Equations 27 and 28 are based on the work of Huang and Nanson (2002) and Huang 
(1996). The parameters a and c are empirically derived and, as indicated in Huang, 
typically have the values a = 4.33 and c = 0.22 pso-" '' , where pso is the median graiil 
diameter at the inlet. The estimated inlet location along the characteristic curve is 
taken to be the point where Equation 27 is satisfied. This relationship was derived for 
river flows, but analogous relationships can reasonably be ^pUed to the channelized 
flows in other types of systems. Ih the alternative, pereons skilled in the art, based on 
the disclosure herem, will recognize o&er inlet conditions could be utilized. For 
example. Equation 27 could be replaced by the assumption that along the 
characteristic curve in the negative-s direction, the deposit thickness will first become 
zero at the inlet.. 

[0048] In tiie sixth step of tiiis embodiment of the method. Figure 5, step 506, the 
extrapolated maxhnum deposit tiiickness obtained in step 504 is checked for 
consistency with the extrapolated inlet flow conditions obtained in step 505. It is 
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assumed that the maximum thickness of a deposit is the height at which the deposit 
would divert the flow that created it In temis of the inlet Froude number. 



Inlet Froude Number: 



the maximum deposit thickness would be. 



=4^ 



[30] 



[0049] Equation 30 derives from the work of Baines (1995). Because Baines' 
work involved two-dimensional air current flows over obstructions, a scaling constant 
q has been added to relate the two-dimensional result to the three-dimensional 
behavior of the hydrodynamic flows of interest. In the present method, scaling 
constant q is determined from numerical modeling to be 0.9. This value of » depends 
on the three-dimensional shape and size of the dq)osit, and the relationship of that 
shape and size to the flow's assumed characteristics, and may be modified for other 
flows. 

[0050] If Equation 30 is not satisfied by the computed maximum deposit 
thickness and inlet flow conditions, then the estimate of flow height in step 501 is 
modified, and steps 502 through 506 are repeated until the condition in Equation 30 is 
satisfied. 

[0051] In the seventh and final step of this embodiment of the method. Figure 5, 
step 507, the inlet conditions calculated when Equation 30 was satisfied are used as 
boundary conditions for solving Equations 1, 2, 5, 6, 7 and 8 to compute the deposit 
size, shape, and internal grain size distribution at each point within the deposit 
Equations 5. 6, 7, and 8 are flow equations used to predict the properties of 
sedimentary bodies. Persons skiUed in the art will recognize many different methods 
to solve Equations 5 through 8 (for example, Bradford & Katopodes, 1999). All 
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methods that can be used to solve these flow equations to detennine the properties of 
flie deposit are intended to be within the scope of this invention. 
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Table 1 



Parameter 


Numeric 
Reference in 
Figures 


Definition of Parameter 


b 


8(17 is the 
iiilet) 


Half-width of the inlet. The inlet is located at x=0 and y 

between H an A 


ho 


26 


Height of the turbid water layer at the inlet, constant across 
the width of the inlet. 


Uo 




Velocity of the turbid water flow out from the mlet. The 
flow is assumed to be in the positive x direction, and the 
T vj.«ur\^x«.jr odouiiicu lu oe Qonsiani laterally across the width 
of the inlet and vertically through the height of the flow. 


Cto 

h(x,y) 
u(x,y) 


27 


Total volume fraction of sediment in the turbid water at the 
inlet. 

Height of the turbid water layer. 

x-velocity component of the turbid water layer 


v(x,y) 
yv(x,y) 
T 

A(x,y,t) 




jr T wyy^xtjr v^umpuiiPxiL Qi luG lurDio. Water laver. 

Magnitude of the flow velocity, w - sqrtfu^ + v^) 

Tmie durmg which the deposit is built before it diverts the 
flow. 

Tune dependent elevation of the subsurface topography 


A(x,y,0) 


24 
30 


■cicvduuu oi tne suDSunace topography before deposition. 
Blevation of the subsurface topography after deposition for 
time period T. 


Zmax 




Maximum tihickness of the deposit. Equal to [A(xni,ynbT) - 
A(xm,ym,0)], where the subscripts indicate the x,y location of 
the maximum thickness of the deposit 


CiCx^) 




Volume fraction of the turbid water layer consisting of 
sediment in the ith grain size bin, averaged vertically through 
the turbid water laver. 


Cj(x^) 




Total volume fraction of the turbid water layer consisting of 
sediment. This is the sum of Ci(x,v) over all grain size bins. 


Gi(x,y) 

Ei(x,y) 
Di(x,y) 
di 




uram size distribution m the deposit, assumed constant 
_ttirough the thickness of the deposit. Sum of G, over all i is 1 
irosion rate for the ith grain size bin 
Deposition rate for the ith grain size bin 


Pso(x,v) 

vsi 

R 




Grain diameter of the i-th grain size bin 

Median gram diameter in the deposit at location x.v 

Settling velocity of grains of diameter di in still water 

Ratio of the density difference between the sediment and 

water to the densitv of water 

Porosity of the deposit. Assumed typically to be 0.4 

Gravitational constant 9.8 m/s^ 

Bntrainment coefficient. 
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f 




Drag coefficient for erosion, typically 0.002. 


CD 




Bottom drag coefficient for flow, typically 0.002. 


V 




Kinematic viscosity of water. Assumed to be 0.01 cm^/s 


s 




Measurement location on the characteristic curve of the flow 


9 




Scaling constant qjplied to predicted maximum deposit 
height, typically 0.9. 


Fr 




Froude number 






Particle Reynolds dimensionless parameter of flow 


C 




Empirical coefficient in the inlet condition 


a 




Empirical coefficient in the inlet condition 


z(x,y) 


36 


Thickness of deposit resulting from deposition 


To 




Ratio of sediment volume fraction at the bottom of the flow 
to Ct. Used in a correction to the deposition rate which 
accounts for vertical variation in the sediment concentration. 
A constant assumed typically to be 2. 
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